Method and Apparatus for Estimating 
The Frequency and/or Phase of a Digital Signal 

BACKGROUND OF THE INVENTION 

[0001] This application claims a priority from German patent application No. 
100 39 666.6, filed on August 14, 2000, and the contents of that application are 
incorporated herein by reference. 

[0002] This invention relates to a method and an apparatus for estimating 
the frequency and/or phase of a digital signal. 

[0003] A method of estimating frequency is disclosed by J. K. Wolf and 3. 
W. Schwartz, "Comparison of Estimators for Frequency Offset," IEEE 
Transactions on Communications, vol. 38, no. 1, January 1990, pages 124-127. 
It is proposed in this article that the phase of a complex digital input signal be 
differentiated and that the differentiated phase be fed to an averaging filter. 
This article shows that the ideal pulse response of the averaging filter is 
parabolic. The parabolic curve of the pulse response of the averaging filter can 
be approximated relatively well by a trapezoidal curve with an increasing range, 
a constant .range and a descending range. The standard deviation of the 
estimation error increases by only about 6% in comparison with using the ideal 
averaging filter with the parabolic pulse response. Therefore, this could be 
termed a sub-optimal frequency estimator. 

[0004] If the filter with the trapezoidal pulse response mentioned in the 
above article were to be implemented directly, a relatively large number of 


multiplications would have to be performed because each sample value within an 
observation interval must be multiplied by a corresponding coefficient. 

[0005] Therefore, it is an object of this invention to provide a method and 
an apparatus for estimating the frequency and/or the phase of a digital input 
signal, which preferably works without logic-gate-intensive multipliers. 

SUMMARY OF THE INVENTION 

[0006] With regard to the method, this object is achieved by: determining 
the phase values of an input signal; summing, or adding up, these input-signal 
phase values over a predetermined summation length N/B which is a 
predetermined fraction 1/B of an observation length of N phase values, to create 
added-up phase values; reducing the sampling rate of the added-up phase 
values by the factor N/B in comparison with a sampling rate of the input-signal 
phase values; delaying the added-up phase values with at least B-l delay 
elements, each of which delays the added-up phase value by one sampling 
period of the reduced sampling rate; adding up the differently-delayed added-up 
phase values to create a resulting pulse response of the frequency so that the 
resulting pulse response of the frequency is constant positive in a first interval, is 
zero in a second interval and is constant negative in a third interval, and/or a 
resulting pulse response of the phase such that the resulting pulse response of 
the phase is constant in at least a middle interval of the observation length and 
is otherwise zero. 

[0007] With regard to an apparatus, this object is achieved by having: a 
phase recording device which determines the phase value of an input signal; a 
first filter, which adds up the input-signal phase values over a predetermined 
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summation length N/B, which is a predetermined fraction 1/B of an observation 
length of N phase values, to form added-up phase values, and reduces the 
sampling rate of the added-up phase values by a factor N/B in comparison with a 
sampling rate of the input-signal phase values; a second filter which delays the 
added-up phase values in a chain of at least B-l delay elements, each of which 
delays the added-up phase values by one sampling period of the reduced 
sampling rate, and adds or subtracts the differently-delayed added-up phase 
values, to create a resulting pulse response of the frequency so that a resulting 
pulse response of the frequency is constant positive in a first interval, is zero in a 
second interval and is constant negative in a third interval, and/or they are 
added to create a resulting pulse response of the phase so that the resulting 
pulse response of the phase is constant in at least a middle interval and is 
otherwise zero. 

[0008] Advantageous refinements and enhancements of this invention are 
also set forth herein. 

[0009] This invention firstly makes use of the recognition that it is more 
advantageous to start with phase values of the input signal rather than phase 
difference values to estimate the frequency and phase. Accordingly, a pulse 
response of an ideal averaging filter varies from a parabolic curve to a curve that 
descends linearly from positive values to negative values. The sub-optimal pulse 
response, which represents a good approximation to the ideal pulse response, 
has then first a positive square-wave range, a range where the pulse response is 
0, and a negative square-wave range. The ideal pulse response for the phase is 
constant over the entire observation length. The component sections with a 
constant course of the pulse response correspond to a summation of the phase 
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values in these component ranges. By a simultaneous reduction of the sampling 
rate 7 the summed phase values are made available at the output of the first filter 
stage timed with, or in cycle with, the fraction of the observation length that 
corresponds to the constant component ranges of the pulse response. In a 
second filter stage, the phase values at this reduced timing are delayed and then 
are added together or subtracted from one another so that the above-described 
characteristic curve of the resulting pulse response is obtained for the estimated 
frequency and/or the estimated phase. 

[OOIO] The first filter stage can be implemented by a series connection of 
an integrator, a sampling rate converter and a following differentiator. In this 
regard, the integrator continuously adds up or integrates the phase values, while 
the differentiator subtracts the starting value at the beginning of the summation 
interval from the final value at the end of the summation interval. 

[0011] In an especially advantageous embodiment of this invention, the 
resulting pulse response for the phase is constant only during a middle interval 
that is 2/3 the total observation length and is otherwise is equal to 0. However, 
the resulting pulse response for the frequency is constant and positive in the first 
third of the observation length, is equal to 0 in the second third of the 
observation length and is constant and negative in the third of the observation 
length. In comparison with a completely constant pulse response over the 
complete total observation length for the phase estimate, this has the advantage 
that the pulse response can be broken down into blocks of one sixth of the total 
observation length and can be added up or subtracted with a suitable 
standardization of these blocks so that the desired resulting pulse response for 
the phase and for the frequency can be created without multiplication, because 
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the result is subject only to a factor of an integral power of the base 2 and thus, 
instead of multiplication, only bit shifting or a shifting of the value of the bits 
must be performed in interpreting the results. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] Two embodiments of this invention are described in greater detail 
below with reference to the drawings. The described and drawn features can be 
used individually or in preferred combinations in other embodiments of the 
invention. The foregoing and other objects, features and advantages of the 
invention will be apparent from the following more particular description of the 
invention, as illustrated in the drawings in which reference characters refer to the 
same parts throughout the different views. The drawings are not necessarily to 
scale, emphasis instead being placed upon illustrating principles of the invention 
in a clear manner. 

[0013] Figure 1 is a block diagram of an apparatus of this invention for 
estimating the frequency and phase of a digital input signal; 

[0014] Figure 2 is a graph of the ideal weighting coefficient w f opt for 
estimating the frequency; 

[0015] Figure 3 is a graph of the sub-optimal weighting coefficient w f so for 
estimating the frequency; 

[0016] Figure 4 is a graph of the optimal weighting coefficient w<p opt for 
estimating the phase; 
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[0017] Figure 5 is a block diagram of a first embodiment of an estimator; 

[0018] Figure 6 is a diagram of a resulting pulse response h f for the 
frequency of the estimator shown in Figure 5; 

[0019] Figure 7 is a diagram of the resulting pulse response h<p for the 
phase of the estimator shown in figure 5; 

[0020] Figure 8 is a block diagram of an embodiment of the apparatus of 
this invention for estimating the frequency and phase in a second embodiment of 
the estimator; 

[0021] Figure 9 is a diagram of a resulting pulse response h f for the 
frequency of the estimator shown in Figure 8; and 

[0022] Figure 10 is a diagram of a resulting pulse response hq, for the phase 
of the estimator shown in Figure 8. 

DETAILED DESCRIPTION OF THE INVENTION 

[0023] Figure 1 shows a basic structure of an apparatus 1 of this invention 
for estimating the frequency f al and/or phase cp a i of a digital input signal x(i). 
The frequency f a i and phase (p a i(T 0 ) of a rotational vector at time point T 0 are to 
be estimated from sample values x(i) of the complex rotational indicator, or 
vector. In this regard, the input signal x(i) may contain additive interference 
n(i), for example in the form of white noise with a Gaussian distribution, which is 
constructively shown added in an equivalent summing unit 2. 
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[0024] In a phase recording, or determining, device 3, phase values C a i(i) 
of the input signal x(i) are determined from the following equation, where they 
are standardized to 2%\ 


C al (i) = angle (x(i))/2^ 


(1) 


[0025] Sample values are obtained by sampling at times tj = i»T a2 and can 
be represented as follows: 


where f a i denotes the frequency of the input signal x(i) to be estimated, and T a2 
is the sampling period of sampling frequency f a2 = 1/T a2 with which the input 
signal x(i) is sampled. N values of an observation length of N phase values C a i(i) 
over an observation time T = N • T a2 are used to estimate the frequency f ai and 
phase (p a i as follows: 


where C a i is a vector (vectors are printed in bold in this patent application) of N 
phase values standardized to 2te. If it is assumed that the interference n(i) is a 
Gaussian-distribution white noise signal, a vector of optimal weighting factors 
w f opt for estimating the frequency f al can be determined by using the maximum 
likelihood theory (see, for example, J. K. Wolf and J. W. Schwartz, 
"Comparison of Estimators for Frequency Offset," loc. at.), where the optimal 
estimated frequency f a i opt can be determined by vector multiplication as follows: 


x(i) = exp ■ (j (2-;r-fai ■ i ■ T a2 + cp 0 ) 


(2) 


C a i = (C al (0),C al (1),..., Cai(N-l)) 


(3) 
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f al ° pt = w f opt ■ c a r 


(4) 


where w f opt is a vector of N weighting coefficients w f opt (i): 

Wf opt = (w f opt (0) , w f opt (1), w f opt (N-l) (5) 


where 


,0Pt (i) = I2.i-6.(N-1) (6) 
(N-1)-(N + 1)-T a2 


[0026] Averaging with the weights w f opt is performed in the frequency 
estimator 5 of the estimator 4 that follows the phase recording device 3. The 
estimated frequency f ai of the input signal x(i) is available at the output of 
frequency estimator 5. 

[0027] Figure 2 shows the weighting coefficients w f opt (i) according to 
equation (6), which illustrates the linear increase in weights with an increase in 
index i. By integrating the function of Figure 2 of the weighting coefficients one 
gets the parabolic characteristic curve of the weighting coefficients given by J. 
K. Wolf and J. E. Schwartz, "Comparison of Estimators for Frequency Offset/' 
loc. at., for the estimate based on the phase difference values. 

[0028] Direct implementation of optimal weighting coefficients given in 
equation (6) and illustrated in Figure 2 has the disadvantage that a great many 
multiplication operations must be performed because the weighting coefficients 
w f opt (i) are different for each sample. In the case of a hardware implementation, 


this leads to a high logic-gate complexity. Therefore, it is proposed that instead 
of the optimal weighting coefficients w f opt (i), sub-optimal, approximated 
weighting coefficients w f so (i) be used, as illustrated in Figure 3. To do so, the 
observation length N is divided into three component ranges, each with a length 
of N/3. In the first third between 0 and N/3-1, the weighting coefficients are 
constant negative; in a second third between N/3 and 2/3N-1, the weighting 
coefficients are constant 0; and in a third one-third between 2/3 N and N-l, the 
weighting coefficients are constant positive. The sub-optimal weighting 
coefficients w f so (i) are then calculated according to the following equation: 

-9/(2 N 2 T a2 ) for 0 < i < N/3 
0 for N/3 < i < 2N/3 (7) 

-9/(2N 2 T a2 ) for2N/3<i<N 

[0029] The standard deviation of the estimation error increases by only 
about 6% in comparison with the use of the optimal weighting coefficient w f opt (i). 
This can therefore be termed a sub-optimal frequency estimator. 

[0030] The phase cp al at the optimal estimation time T 0 opt = (N-l) • T a2 /2 
can also be estimated with a phase estimator 6 of the estimator 4 by multiplying 
the vector C a i times the weighting vector w f opt as follows: 

<PaiO"o opt ) = w, 0 '*-C.i T (8)] 

where 

w/"* = ( W(p opt (0) , w, opt (1) , w, opt (N-l)) (9)] 


w f so (i) = 
w f so (i) = 
w f so (i) = 
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[0031] For the weights w<p opt (i) the following holds: 


w^ 1 * (i) = 2*/N (10)] 

[0032] The optimal weighting coefficients W(p opt (i) for estimation of the 
phase cp a i are thus constant over the entire observation length N. An 
approximation by means of sub-optimal weighting coefficients is not necessary 
for the phase because of this detail. The weighting coefficients Wq, opt (i) are 
shown in Figure 4 as a function of the index i. 

[0033] Figure 5 shows an implementation of the estimator 4 according to a 
first embodiment of this invention. 

[0034] The estimator 4 is roughly divided into a first filter 7 and a following 
second filter 8. The first filter 7 is a cic (cascaded integrator comb) filter of the 
first order. The first filter 7 includes an integrator 38 with an adder 9 and a 
delay element 10, a differentiator 13 with a delay element 11 and a subtracter 
12, as well as a sample-rate converter 14 (down-sampler) arranged between the 
integrator 38 and the differentiator 13 or difference former. The sample-rate 
converter 14 reduces the sampling rate f a2 of the phase values C a i at the input 
of the first filter 7 by the factor N/3 in this embodiment, where N stands for the 
observation length, that is, the number of observed phase values C a i(i) during 
the observation time N • T a2 . The integrator 10 continuously adds up various 
phase values C a i(i) supplied to it. The sampling-rate converter 14, in 
combination with the differentiator 13, limits the summation length to N/3, 
because the differentiator 13 subtracts the starting value at the beginning of the 
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summation from the end value of the summation over N/3 phase values C a i(i). 
Therefore, after N/3 sampling periods T a2 of the phase values C a i(i), a sum value 
which is referred to below as the added-up phase value S a i(i), is obtained at the 
output of the first filter 7 and represents the sum of the previous phase values 
Cai(i). 

[0035] These added-up phase values S a i(i) are fed to the second filter 8. 
In this embodiment, the second filter 8 includes two delay elements 15 and 16 
which delay the added-up phase values S a i(i) in a timing cycle (3/N) • f a2 by one 
sampling period N • T a2 /3 each. In this regard, the added-up phase value S a i(i) 
is fed to the "+" input of a subtracter 18, whereas the added-up phase values 
S a i(i-2) which has been delayed by two sampling periods is fed to the "-" input of 
subtracter 8, so that the phase value S a i (i) is delayed by the delay time 2 • N • 
T a2 /3 and appears inverted at the output of the subtracter 18. 

[0036] The cic filter 7 creates a square-wave pulse response with a 
constant positive coefficient over the length N/3. If the inverted pulse response 
is added to this with a delay of 2 • N/3, the resulting pulse response h f shown in 
Figure 6 with an amplitude of 1 is obtained. By multiplying by the factor 9/(2 • 
N 2 • T a2 ) in a multiplier 19, the resulting pulse response h f shown in Figure 6 is 
obtained for the frequency over the observation length N. Since the estimated 
value for the frequency is of interest only in the periodicity of the observation 
length N, in another sampling-rate converter 37 the sampling rate is reduced by 
a factor of 3 again, i.e., on the whole by a factor of N in comparison with 
sampling rate f a2 . By mirroring the pulse response h f shown in Figure 6 on the 
ordinate which is shifted by N/2, the function of the sub-optimal weighting 
coefficient w f so for the frequency as shown in Figure 3 is obtained exactly. This 
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mirroring on the offset ordinate corresponds to the different definition of the 
vector multiplication on which equation (4) is based in comparison with the 
definition of convolution or folding of the input sequence C al (i) with the resulting 
pulse response h f . As is known, the output function of a filter is obtained by 
folding the input function of the filter with the pulse response of the filter, which 
is equivalent to vector multiplication of the input sequence by the mirrored 
coefficients of the pulse response. 

[0037] The second filter 8 has two adders 20 and 21 for creating the 
resulting pulse response h^ for estimation of the phase cp a i. The added-up phase 
value S a i(i) at the input of the first delay element 15, which adds up the added- 
up phase value S al (i-1), which is delayed by (N/3) • T a2 at the output of the first 
delay element 15, and the added-up phase values S ai (i-2) which is shifted by (2 
• N/3) • T a2 . As mentioned above, the cic filter 7 creates a positive constant 
(square-wave) partial pulse response of the length N/3. By adding three such 
partial pulse responses of the length N/3, this yields on the whole positive 
constant pulse response over the total observation length N. By multiplication by 
the factor 2rc/N in the multiplier 22 and reducing the sampling rate by the 
remaining factor 3 in the sampling-rate converter 23, the resulting pulse 
response h Q for the phase is obtained as shown in Figure 7. 

[0038] However, one disadvantage in the embodiment illustrated in Figure 
5 is that the factors 9/(2 • N 2 • T a2 ) and 2te/N of multipliers 19 and 22 are not an 
integral power of the base 2. Therefore, true multiplication is necessary in 
multipliers 19 and 22, not just a shift in the value of the bits of estimated 
frequency f a i which is represented in binary notation and estimated phase cp ai , 
necessitating a relatively great logic-gate complexity in implementing the 


12 


circuitry. However, it should be emphasized that the gate complexity in the 
embodiment of Figure 5 is clearly reduced in comparison with use of optimal 
weighting factors h f opt shown in Figure 2, because it is only necessary to perform 
multiplication after N/3 sampling periods T a2 . 

[0039] In the embodiment illustrated in Figure 8, these remaining 
multiplications are also eliminated or replaced by bit shifting. In the embodiment 
illustrated in Figure 8, elements with the same reference notations corresponding 
to those in the embodiment illustrated in Figure 5 are the same, so it is not 
necessary to repeat the description. 

[0040] First, the phase cp a i is standardized to In. Furthermore, the 
frequency estimate value f ai is standardized to the inverse of the sampling period 
1/T a2 . In this embodiment, a counter 24 which counts the frequency f ai of the 
input signal x(i) continuously is used to create the phase values C a i(i) from the 
input sequence x(i). The counter status of counter 24 is sampled or read out at 
times t| = I»T a2 , which is represented by switch 25. The phase values C a i(i) thus 
created are fed to the cic filter 7. Even better results can be obtained with a 
dual-flank counter. It should be emphasized that counter 24 is only one of the 
many possibilities of obtaining phase values C a i(i). 

[0041] In the embodiment illustrated in Figure 8, the sampling-rate 
converter 14 reduces the sampling rate f a2 = 1/T a2 by the factor N/6, so that cic 
filter 7 adds up the phase values C a i(i) only over the fraction 1/6 of the 
observation length N. The added-up phase values S a i(i) are fed to a chain of five 
delay elements 26, 27, 28, 29 and 30 in the second filter 8. Since the cic filter 7 
creates only a constant, positive pulse response of the length N/6, first two 
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positive partial pulse responses immediately in succession must be combined to 
create the resulting pulse response h f for the frequency shown in Figure 9 and 
then, delayed by N/3, two inverted partial pulse responses must be combined 
again in immediate succession. This can be achieved because phase values 
S a i(i) at the input of the first delay element 26 and phase values S a i(i-1) at the 
output of the first delay element 26 are fed to an adder 31 to create the positive 
component of the resulting pulse response h f between 0 and N/3-1. The outputs 
of the second delay element 27 and third delay element 28 do not enter directly 
into the formation of the resulting pulse response h f , so that the resulting pulse 
response h f is identical to 0 between N/3-1 and 2N/3. The added-up phase 
values S a i(i-4) at the output of the fourth delay element 29 and S a i(i-5) at the 
output of the fifth delay element 30 are respectively fed to the minus "-" inputs 
of subtracters 32 and 33, with the output of the adder 31 being fed to the plus 
"+" input of subtracter (which could also be termed a negative adder) 32 and 
the output of subtracter 32 being fed to the plus "+" input of subtracter 33. 

[0042] Since the value was standardized to 1/T a2 the multiplier is now 9/(2 
• N 2 ). In the sampling-rate converter 37, the sampling rate must be reduced by 
the remaining factor 6. 

[0043] In the embodiment illustrated in Figure 8, a reduced effective 
observation length of 2/3 N is used to estimate the phase <p a i. To do so, the 
resulting pulse response l\p for the phase is set at 0 in the interval between 0 
and N/6-1 and between 5/6N and N-l by not including the added-up phase 
values Sai(i) at the input of the first delay element 26 and S a i(i-5) at the output 
of the fifth delay element 30 in the addition to create the resulting pulse 
response h^. The added-up phase values, S ai (i-1) at the output of the first delay 
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element 26, S al (i-2) at the output of the second delay element 27, S ai (i-3) at the 
output of the third delay element 28, and S a i(i-4) at the output of the fourth 
delay element 29, are added in the adders 34, 35 and 36. Because of the 
standardization to 2n and the reduction of the effective observation length from 
N to 2/3 N, this yields a correction factor of 3/2 N. In the sampling-rate 
converter 23, the sampling rate for the phase must still be reduced by the 
remaining factor 6. 

[0044] The resulting pulse response h^ for estimating the phase cp ai is 
shown in Figure 10. 

[0045] The main advantage of the embodiment illustrated in Figure 8 in 
comparison with the embodiment illustrated in Figure 5 is that the two correction 
factors 9/(2N 2 ) and 3/(2N) can be created by bit shifting. This is true because if 
N = 3»2 n , where n is a whole positive number, then 9/(2 • N 2 ) = 2" 2n l . This is 
division by a power of 2 which, when manifested as hardware, can be 
represented as a displacement (shifting) operation. The same thing then also 
applies for the factor 3/(2N), because with N = 3«2 n it holds that 3/(2»N) = 2""" 1 . 
This is also division by a power of 2, which can be represented as a bit shift 
operation in a hardware implementation. It is therefore possible to eliminate the 
multipliers 19 and 22, which are shown in Figure 8 only for the sake of 
illustration and replace them with a reinterpretation of the value of the bits at 
the output of the subtractor 33 or adder 36. 

[0046] This invention is not limited to the embodiments shown herein. 
Instead of a summation length of N/3 or N/6 in the cic filters 7, a fraction 1/B of 
the observation length N may also be used in general as summation length N/B, 
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wherein then the sampling rate f a2 must be reduced by the factor N/B in the 
sampling-rate converter and by factor B in the sampling-rate converters 20 and 
23. To eliminate multipliers 19 and 22, it is advantageous if this fraction is 1/B = 
1/(3 • n), where n is a whole positive number. 


16 


